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This  paper  is  a  continuation  of  a  preliminary  study  which  investigated  the 
magnetic  field  radiated  from  a  charge  bunch  traveling  over  a  finite  path.  One  of  the 
recommendations  of  the  preliminary  study  was  to  investigate  how  the  magnitude  and 
shape  of  the  radiated  field  depends  on  the  location  of  the  observer  relative  to  certain 
time  boundaries  defined  by  relationships  between  the  arrival  times  of  different  parts  of 
a  pulse.  This  recommendation  was  followed  by  studying  the  signals  radiated  from 
charge  bunches  with  two  different  ratios  of  rise  time  to  path  length.  Within  the 
Cerenkov  cone  and  near  a  time  boundary,  the  observed  pulses  are  shaped  like  spikes. 
The  time  boundaries  were  found  to  be  hyperbolas  in  the  plane  defined  by  the  direction 
of  travel  of  the  charge  bunch  and  the  direction  of  propagation  of  radiation. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  topic  of  time  developement  for  Cerenkov  radiation  was  introduced  by 
Buskirk  and  Neighbours  in  [Ref.  1]  where  the  concept  was  applied  to  an  electron  beam 
with  a  ramp-front  function.  Later  in  [Ref.  2],  which  was  a  preliminary  study,  magnetic 
fields  were  calculated  and  plotted  with  respect  to  time  for  various  locations  relative  to 
the  beam.  The  values  of  p'  and  the  index  of  refraction  n,  were  chosen  in  order  to  study 
the  relationship  between  the  time  regions  relative  to  the  beam  and  the  shape  of  the 
radiated  magnetic  field  curves.  The  value  P'  is  equal  to  v/c,  where  c  is  the  velocity  of 
light  in  the  medium.  From  this  preliminary  study,  one  of  the  recommendations  was  to 
investigate  how  the  magnitude  and  shape  of  the  radiated  field  depends  on  the  location 
of  the  observer  relative  to  certain  time  boundaries  defined  by  relationships  between  the 
arrival  times  of  different  parts  of  a  pulse. 

The  B  field  calculations  in  [Ref.  2:  pp.  70-79]  are  theoretical  and  have  assumed  a 
linear  ramp-front  charge  profile,  a  one  dimensional  current  distribution,  and  a  rigid 
charge  shape.  In  other  words,  the  radiation  is  produced  by  a  bunch  of  electrons  of 
negligible  transverse  dimension  with  a  longitudinal  distribution  which  does  not  change 
as  the  bunch  travels  along  the  z  axis.  The  radiated  power  is  proportional  to  the  square 
of  the  charge  per  unit  length  and  therefore  to  the  square  of  the  beam  current.  The 
expressions  for  B  and  E  in  [Ref.  1:  p.  3751]  are  valid  and  the  evaluation  is  therefore 
coherent  only  insofar  as  the  bunch  is  not  distorted  either  by  the  reaction  of  the 
radiation  or  by  the  instabilities  associated  with  very  high  current  beams,  [Ref.  1:  p. 
3753].  So  far  none  of  the  calculations  have  been  compared  to  experimental  results  due 
to  the  shortage  of  suitable  data. 

Present  technology  is  such  that  this  structure  is  not  observable  in  the  Cerenkov 


radiation  from  S  or 


band  linacs  because  of  their  relatively  high  fundamental 


frequency.  However,  other  accelerators  with  their  longer 


electron  bunch  structure  and 


higher  currents,  should  produce  observable  Cerenkov  signals  in  air  for  energies  greater 
than  about  25  MeV.  [Ref.  1:  p.  3753] 


B.  OBJECTIVES 
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The  main  objective  of  this  work  is  to  investigate  the  magnitudes,  shapes,  and 
pulse  widths  of  the  radiated  fields  as  a  function  of  position  in  the  radiation  plane. 
More  specifically,  this  thesis  follows  the  recommendation  to  investigate  the  influence  of 
the  time  boundaries  in  this  plane.  These  time  boundaries  are  known  to  depend  on  the 
arrival  times  of  signals  from  different  parts  of  the  charge  distribution. 

The  secondary  objectives  of  this  thesis  are  a  result  of  the  main  objective.  In  order 
to  calculate  the  B  field  as  a  function  of  time,  the  FORTRAN  program  of  [Ref.  2:  pp. 
70-79]  needs  to  be  written  in  a  structured  format  and  transtered  to  the  NFS  main  frame 
computer.  Appendix  A  contains  this  program  called  "CERENKV"  which  is  written 
using  the  WF77  version  of  FORTRAN.  This  program  evaluates  the  integral  which 
yields  the  magnitude  of  the  Cerenkov  pulse  at  a  given  time  after  an  electron  beam  is 
discharged  into  a  medium  from  an  accelerator.  The  last  objective  provides  yet  another 
tool  needed  for  studying  the  characteristics  of  the  radiated  magnetic  field.  This  is  a 
FORTRAN  program  that  will  generate  the  data  for  the  time  boundaries  which  plot  the 
S/L  verses  the  Z/L  (S-Z  plane)  graph.  In  the  next  chapter,  S  is  defined  as  the  radius 
vector  and  Z  is  the  horizontal  distance  measured  from  the  beam  discharge  point  along 
the  beam  path  length  L.  The  program  is  found  in  Appendix  B  and  is  available  on  the 
NPS  main  frame  under  the  filename  "TIME”. 
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II.  THEORY  AND  EQUATIONS 

a.  Radiation  fields 

The  potentials  from  a  moving  charge  distribution  were  determined  in  [Ref.  1:  p. 
3750).  The  fields  (in  cgs  units)  from  these  potentials  are  given  by  Equations  2.1  and 


B  =  V  x  A 


E  =  -  V<D  -  (1  co)(d\  dt) 


The  magnetic  field  B,  may  be  calculated  from  Equation  2.1  and  since  A  has  only  a  z 
component,  this  approach  is  easier  than  finding  E.  A  similar  derivation  can  be  made  in 
order  to  find  the  magnitude  of  the  E  field.  It  is  also  true  that,  in  the  Cerenkov  case, 
E  B  =  c  cQ,  which,  for  plane  waves,  is  the  usual  relation  between  the  electric  and  and 
magnetic  fields.  [Ref.  1:  p.  3752]. 

In  [Ref.  1:  p.  3750],  it  was  assumed  that  a  charge  density  function  pv  and  a 
current  density  j  =  pv  cq  moved  with  velocity  v  in  the  plus  z  direction.  The  charge  and 
current  are  assumed  to  be  concentrated  along  the  z  axis  such  that 

p0(r,t)  =  p(z.t)6(x)6(v).  (2.3) 

The  charge  is  assumed  to  move  with  no  change  in  shape  so  that  the  z  and  t 
dependence  of  the  charge  is  p(z,t)  =  po(z-vt).  The  usual  charge  and  current  densities 
are  represented  by  pir  and  jv.  while  p  and  po  are  charge  per  unit  length.  The  velocity  of 

' : 5 h t  :  •  :  and  :  w  *he  medium  end  free  mace  -esrectivelv. 

The  potential  ire  found  n-  ui.-uug  *he  usual  retarded  solutions  to  'he  wave 

equations;  which  become  under  the  assumption  of  a  line  distribution  of  charge 
(Equation  2.3), 


Afr.tl  -  (vco)f(I  R)p(r',t)dz'. 
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where  R  =  r  -  r'  and  t\  the  retarded  time,  is  given  as  t'  =  t  -  |r  -  r'|/c.  Introducing  a 
new  variable,  u(z')  =  z'  -  vt'  into  Equation  2.4  gives 

A(r,t)  =  (v/co)J(l/R)po(u)dz'.  (2.5) 

Since  the  charge  is  confined  to  the  z'  axis,  the  new  variable  u(z')  can  be  written  more 
explicitly  as 


u(z  )  =  z'  -  vt  +  (v/c)[x2  +  y2  +  (z  -  z')2]'’2. 


(2.6) 


The  potential  A,  has  only  a  z  component,  and  the  B  field  has  only  x  and  y 
components,  Bx  =  dAJdy  and  By  =  -  dAJdx.  Carrying  out  the  differentiation  for  the 
x  component  gives 

Bx  =  v/cJ(dR'I/5y)po(u)dz'  +  v/cjR'1(ap0(u)/ay)dz'.  (2.7) 

For  radiation,  the  first  integral,  falling  off  as  R'2  at  large  distances,  will  be  neglected 
and  only  the  second  term  is  considered  further.  From  Equation  2.6,  it  is  seen  that  u  is  a 
function  of  x  and  y  so  that  the  second  integral  can  be  written 

Bx  =  v2/ccJ(y/R2)p'o(u)dz',  (2.S) 


where  p  Q(u)  is  the  derivative  of  po  with  respect  to  its  argument  u.  The  corresponding 
expression  for  By  has  y  replaced  by  -x.  These  two  components  combined  give  the  total 
magnetic  radiation  field  B.  In  the  cylindrical  coordinates  (s,0,z),  where  s  is  the  radius 
vector,  s  =  (x2  +  y2)'/2,  B  is  tangential  (i.e.,  in  the  0  direction)  with  a  magnitude 
given  by 


3  =  .'v2,oc  )f{s/RV  u)dz*. 

•>'j '  r  >' 


Writing  Equation  2.9  with  the  finite  limits  z'{  and  z'f  with  v2/ccq  =  nf2  and  assuming 
p'(u)  =  pm  =  constant,  the  new  expression  for  B  is  written  as 


Here,  s  -  (x:  +-  >•*)'*  and  R:  -  s*  +  (z  -  z  )"  or  R;  *  s:  -  w*.  where  w  -  (z  -  z). 

Prior  to  evaluating  Equation  2.10,  it  is  necessary  to  consider  the  dependence  of  u 
on  z  as  well  as  p '  on  u.  As  previously  stated,  the  charge  density  function,  p(z.t) 
moves  with  velocity  v  in  the  plus  z  direction  without  changing  shape.  This  charge 
density  motion  makes  up  tr.e  beam  current  profile.  The  upper  drawing  ::t  Figure  2.1 
illustrates  the  charge  density,  p(z.t).  The  nse  in  the  charge  density  is  approximated  to 
be  'linear  until  it  reaches  a  maximum  constant  value  and  then  declines  linearly  to  zero. 
The  time  it  takes  for  the  charge  to  reach  this  plateau  is  referred  to  as  the  rise  time. 
Multiplying  the  rise  time  bv  v  Yields  the  value  Au  =  (u.  -  u,)  which  is  used  later  in 
this  paper.  The  lower  drawing  is  the  derivative  of  the  charge  density  function  which  is 
the  constant,  p  m  in  Equation  2.10. 


Since  u(z')  and  P’0iu)  are  coupled  together  by  a  common  axis  u,  the>  were  closely 
compared  side  by  side  m  [Ref.  1:  p.  3'r51J.  Figure  2.2  illustrates  this  comparison  and 
aids  as  a  supplement  for  the  discu^on  which  follows.  The  first  two  terms  in  Equation 
2.6  are  a  straight  line  with  unit  slope  and  an  intercept  which  changes  with  time,  while 


the  third  term  is  a  hyperbola  opening  in  the  +  u  direction  with  asymptotic  slopes  of 

iv  c.  The  sum  of  these  two  curves  is  u(z).  The  limiting  slopes  are  1  ±  v.'c.  Here,  1  + 

v  c  is  approximately  2.  while  1  -  v/c  is  close  to  zero.  In  the  Cerenkov  case,  1  -  v.  c  <  0 
and  all  the  curves  for  u(z')  for  any  time  t,  have  a  minimim.  With  v  >  c,  the  resultant 
ui,z  )  is  a  curve  whose  ends  both  point  upward  as  shown.  As  time  increases,  the  entire 
curve  will  translate  downward  to  smaller  u  values  as  a  result  of  the  negative  term  (-vt) 
in  Equation  2.6.  For  the  purposes  of  this  paper,  the  Uj  value  will  be  the  zero  starting 
point  and  the  u,  value  will  be  a  negative  value  when  entering  data  in  the  program  of 
Appendix  A. 


Figure  2.2  The  Function  u(z')  and  p'  (u). 


Only  changing  currents  (those  with  a  nonzero  p'Q)  will  contribute  to  the  magnetic 
:n  .ram  Equation  2.1°.  Again,  vnce  me  end  ;f  the  current  puise  increases 
..near.)  upward,  tne  derivative  p0(,uj  wiil  be  a  constant  valued  squared  puise  of 
magnitude  p  m  and  is  also  shown  in  Figure  2.2.  The  corresponding  negative  p  0(u) 
pulse  occuring  at  the  tail  of  a  current  pulse  is  not  shown  and  its  effect  is  considered 
separately.  At  a  time  t.  where  the  u(z')  curve  is  tangent  to  the  upper  portion  of  the 
p  0(u)  pulse,  the  B  pulse  begins.  The  value  of  the  integral  in  Equation  2.10  increases  as 


u(z')  continues  its  constant  downward  motion  with  increasing  time  until  it  becomes 
tangent  with  the  lower  part  of  the  p  0(u)  pulse.  At  this  time  the  nonzero  part  of  the 
integral  has  the  largest  extent.  At  later  times,  the  integral  will  break  up  into  two 
regions  of  the  z'  axis  and  if  p'0(u)  is  constant,  the  value  of  the  integral  decreases  with 
increasing  time  because  the  extent  of  the  integral  in  the  two  regions  continues  to 
decrease  as  a  result  of  the  upward  turn  of  u(z').  The  preceding  discussion  applies  to  z' 
ranging  from  -°o  to  +  °o  with  Uj  and  u2  serving  as  the  upper  and  lower  bounds  in  the 
integration  process,  [Ref.  1:  p.  3751]. 

The  time  structure  (shape)  of  the  resulting  B  pulse  can  be  determined  by 
integrating  Equation  2.10  directly,  [Ref.  1:  p.  3751].  This  is  where  the  limits  of 
integration  must  be  defined  explicitly.  The  limits  in  the  preceding  paragraph  are  found 
with  relative  ease  if  the  bounds  of  Uj  and  u2  extend  indefinitely  along  the  z'  axis  and 
are  the  only  bounds  to  consider.  However,  this  is  not  the  case  considered  in  this  paper. 
Along  with  the  bounds  considered  above,  there  is  a  definite  set  of  boundaries  that  exist 
as  a  result  of  using  a  finite  beam  path  length.  Figure  2.3  shows  the  three  examples  of 
where  the  finite  path  may  fall  in  the  u-z  plane.  These  three  cases  will  be  discussed  in 
more  detail  in  the  next  subsection  and  the  limits  of  integration  will  be  found  for  each 
of  the  three  cases.  Before  preceding  however,  it  is  a  good  time  to  evaluate  the  integral 
in  Equation  2.10. 

Evaluating  the  integral  in  Equation  2.10  is  easy  as  long  as  p'o  is  flat  (or 
constant).  This  assumption  was  made  earlier  in  Figure  2.1.  Using  the  expression  w  = 
(z  -  z')  in  Equation  2.10  and  integrating  with  respect  to  w,  the  final  expression  for  B  is 
found  in  Equation  2.11. 

B  =  np2p  m[arctan(Wj/s)  -  arctan(w2/s)]  (2.11) 

With  the  expression  for  B  now  available,  the  next  step  is  to  find  the  limits  of 
integration  by  finding  the  values  of  Wj  and  w2.  These  limits  were  carefully  derived  in 

rRef.  2:  pp.  15-31].  However,  the  following  sections  will  display  the  key  figures  and 
equations  in  order  to  better  familiarize  future  thesis  students  with  the  terms  and 

physical  concepts  of  Cerenkov  radiation. 


Figure  2.3  Function  u(z')  and  Three  Finite  Paths. 


B.  FINITE  BEAM  PATHS 

1.  Possible  Observing  Positions  P(z,s) 

When  the  path  of  the  beam  is  infinite,  the  limits  of  integration  are  determined 
by  the  intersection  of  the  u(z')  curve  with  the  bounding  values  ut  and  m.  In  the  case  of 
a  path  of  finite  length,  the  integration  limits  are  determined  by  these  intersections,  and 
in  addition,  the  intersection  of  the  u(z')  curve  with  the  limits  of  the  path  length. 

Again,  Figure  2.3  shows  the  three  possible  finite  length  beam  paths,  each  with 
a  different  orientation  with  respect  to  the  minimum  in  u(z').  Path  "a"  is  to  the  left  of 
the  minimum,  path  "b"  contains  the  minimum,  and  path  "c"  is  to  the  right  of  the 
minimum.  These  three  situations  are  labeled  "path  to  the  left",  "path  centered",  and 


"oath  to  the  ri>zht"  respectiveiv. 


nermcr 


three  situations  also  correspond  to 


three  positions  of  interest  relative  to  the  Cerenkov  angle  0C,  where  0C  =  arccos(l  P'). 
This  relation  for  0„  is  called  the  "Cerenkov  relation".  It  is  found  bv  usins  Huveens 
construction  found  in  (Ref.  2:  p.  7].  along  with  basic  trigonometry.  The  actual  relation 
given  was  cos  0  =  I ,  Pn,  where  n  is  the  index  of  refraction  and  B  is  the  dimensionless 


constant  v,  c.  The  latter  is  a  function  of  the  energy  level  and  is  found  by  using  the 
expression  for  the  Lorentz  contraction  factor.  Finally,  the  product  nP  is  defined  as  P'. 

The  three  positions  of  interest  are  relative  to  the  Cerenkov  angle  and 
subsequently  the  beam  path.  The  path  to  the  right  means  the  observer  is  outside  the 
Cerenkov  cone.  The  path  to  the  left  means  the  observer  is  within  the  Cerenkov  cone, 
and  centered  means  the  observer  is  on  the  cone.  These  three  possible  situations  are 
shown  in  Figures  2.4  to  2.6,  which  were  extracted  from  [Ref.  2:  pp.  16,17],  With  beam 
length  L,  and  the  position  of  the  observer  P(z,s),  Figures  2.4,  2.5,  and  2.6  illustrate  the 
path  to  the  right,  to  the  left,  and  centered  respectively. 


Figure  2.4  Path  to  the  Right. 

For  the  path  to  the  right.  0j  must  be  greater  than  0.;  for  the  path  to  the  left.  0,  must 
be  less  than  0c;  for  the  path  centered  about  the  minimum,  0j  must  be  less  than  0c  and 

B,  must  be  greater  man  or  euuai  :o  0 

The  integration  limits  are  determines  by  the  intersections  of  the  u^z  )  curve 
with  the  boundaries  of  the  rectangle  formed  by  the  beam  path  and  the  length  of  the 
pulse  rise.  Changes  in  the  limits  occur  when  the  u(z')  curve  meets  the  comers,  or  when 
it  becomes  tangent  to  the  upper  and  lower  bounds.  The  latter  happens  in  only  the 
"centered"  case.  Labeling  the  corners  of  the  rectangle  as  a,  b,  c,  and  d,  the  boundary 
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Figure  2.5  Path  to  the  Left. 

times  are  found  by  substituting  the  appropriate  boundary  coordinates  into  Equation 
2.6.  Different  situations  arise  depending  upon  whether  diagonally  opposite  or  adjacent 
comers  are  simultaneously  intersected  by  the  u(z')  curve. 

Figure  2.7  shows  successive  intersections  of  the  descending  u(z')  curve  for  the 
path  to  the  right  of  the  minimum.  In  this  case,  t  is  aiwavs  the  earliest  time  and  t.,  is 
the  latest  time.  The  order  of  the  two  times  depends  on  the  relative  size  of  the  beam 
path  length  and  the  length  of  the  pulse  rise.  Figure  2.7  is  drawn  for  t  >  t. . 

The  boundary  times  ta,  tb,  tc,  and  td  are  given  in  Equations  2.12  to  2.15  and  are  found 
in  [Ref.  2:  p.  15]. 


ta  =  [P'(s2  +  z2)Vz  -  u,]  v 

ib  =  (P'(s2  +  z2yl  *  U21  v 


*  [L  +  P'(s2  +  (z  -  L)2),/j  -  u2].v 


(2.121 

(2.13) 


(2.14) 


ft 


Figure  2.6  Path  Centered  About  the  Minimum. 


-  [L  +  P'(s2  +  (z  -  L)2)’  j  -  u2]'v 


(2.15) 
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2.  Limits  of  Integration 

Using  the  critical  times  from  the  previous  subsection,  the  limits  of  integration 
for  Equation  2.11  can  now  be  found.  Since  this  thesis  uses  the  results  of  [Ref.  2:  pp. 
15-31],  the  derivations  are  omitted.  However,  for  the  sake  of  clarity,  Figures  2.S  thru 
2.10  from  (Ref.  2:  pp.  19-21],  are  shown  as  a  review  to  the  way  the  magnetic  fields  are 
calculated  in  Appendix  A. 

As  previously  stated,  Figure  2.7  illustrates  the  relationship  between  the 
function  u(z’)  and  the  finite  path  which  falls  to  the  right  of  the  minimum.  More 
specifically.  Figures  2.8.  2.9.  and  2.10  shows  how  the  limits  of  integration  are  derived  as 
some  arntrarv  time.  :  oasses  tnrcueh  the  critical  times  t. .  ,  t..  and  t ..  In  i  Ref.  2:  ro. 
19-20],  the  expressions  for  z  -  and  z'f  were  derived  with  the  new  variables  Aj  =  Uj  +  vt 
and  A,  =  u2  +  vt.  Equations  2.16  and  2.17  arc  the  expressions  for  z';  and  z'f 
respectively. 


{(P  2z  -  AO  +  P  [(z  -  A,)2  -  s2(P'2  -  I)],/j}/(P'2  -  1) 


(2.16) 


.’.SV 


Fieure  2.7  Path  to  the  Right:  t  >  t.  . 

*■  ~  C  D 


2/  =  {(P'-Z  -  A,)  +  P'[(Z  -  A,)2  -  S2(p'2  -  I )] 1/2 }  (P'2  -  1) 


(2.17) 


In  addition  to  tc  >  tb,  there  are  two  other  cases  to  consider,  namely  tb  >  tc  and  tb  = 
t..  Figures  2.11  and  2.12  illustrate  these  cases. 

The  path  to  left  case  is  similar  to  the  path  to  the  right  case,  but  now  the  three 
cases  are  ta  >  td<  tfl  <  td,  and  ta  =  td>  These  cases  are  shown  in  Figures  2.13,  2.14, 
and  2.15.  Solving  for  the  limits  of  integration  is  similar  to  the  path  to  the  right,  but 
now  the  negative  part  is  used  due  to  the  position  relative  to  the  minimum  of  u(z'). 
[Ref.  2:  p.  25].  Equations  2.18  and  2.19  give  the  "path  to  the  left"  version  of  Equations 
2.16  and  2.17. 


zi  ~  d  P  ”z  *  a i .*  ■  ?  IvZ  •  Aj--  *  s-p  ■  -  i  )j  *;  vp  ■  - 


(-■is; 


2f  =  {(P'2z  -  A,)  -  p-[(z  -  A,)2  -  S2(p'2  -  1 )] 1 2 } . (P'2  -  1) 


(2-19) 
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Figure  2.8  tc  >  tb  :  ta  <  t  <  tfe. 

When  centered  about  the  minimum,  Figure  2.16  introduces  the  new  variables 
tj,  t„  z'c.  and  t3.  Equations  2.20,  2.21,  and  2.22  are  the  expressions  for  these  variables 
and  are  found  in  [Ref.  2:  pp.  2S.74].  In  Equation  2.22.  f  equals  tan(0c)  and  the  variable 
t,  is  the  larger  of  tfa  or  tj. 


=  u: 


+  P  [s2  +  (z  -  z'c)2]Yl  -  u j } / v 


(2.20) 


t,  - 


{z* 


P  (s~  +  (z  -  z'c)2]t/x  •  u2)'v 


(2.21) 


In  the  centered  case,  the  critical  times  to  integrate  between  are  t}.  t, .  and  t?. 
The  arbitrary  time  t,  is  integrated  between  tj  and  t;  and  then  t2  and  t3.  Figures  2.17 
and  2.18  illustrate  this  integration  process.  The  new  limits  from  Figure  2.17  are  given 
in  Equations  2.23  and  2.24. 
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(2.23) 


Fiaurc  2.9  t  >  t.  :  t.<  t  <  t  . 

-  cob  c 

z'  +  -  {(P  2Z  -  Aj)  +  P'[(z  -  A^2  -  S2(p':  -  1)]’*};(P'2  -  1) 
z*.  =  {(P'2z  -  Aj)  -  P'((z  -  Aj)2  -  s2(P'2  -  1)]'  :}  (p-2  -  1)  (2.24) 


From  Figure  2.18,  the  new  variables  that  must  be  considered  here  are  z'JUj).  z'  +  (Uj), 
z'_(u->),  and  z  +  (u2).  A  •  stated  in  (Ref.  2:  p.  31],  z'_(u,)  is  found  using  Equation  2.24 
for  the  first  integral  and  z'f  will  always  be  z'_(u2).  where  z'Ju2)  is  given  by  Equation 
2.25. 


z'.(u2)  =  ((p  ‘z  -  A,)  -  P’[(z  -  A2)>  -  s2(p2  -  1 )] " 2 } . (P'2  -  1) 


(2.25) 


•,Y’; 


...  •  !  .2  wt  «.  s.  .  .  i  l cl .  .  ow.Wi 


>.  otherwise  z  will  euual  zero. 


Similarly,  the  limits  for  the  second  integral,  z'  +  (Uj)  and  z'^.(u,)  will  be  the 
same  as  the  above  equations  in  which  the  second  term  is  positive.  In  this  case.  z';  will 
always  be  equal  to  z'  +  (u2>.  whereas  z'f  will  equal  z’^(Uj)  if  the  latter  is  less  than  L, 
otherwise  z  f  will  equal  L. 


Figure  2.10  t  >  t.  :  t  <  t  <  t,. 
w  c  o  c  a 

3.  Normalized  Coordinates 

As  a  prelude  to  the  summary  of  results,  it  was  decided  that  working  in 
normalized  time  would  make  the  graphs  produced  by  Appendix  A  easier  to  study.  A 
good  base  time  for  normalization  is  the  time  it  takes  for  an  electron  bunch  to  travel  the 
distance  L.  To  this  end,  dividing  the  boundary  times  (Equations  2.12  to  2.15)  by  x  = 
L'v  gives  the  new  normalized  boundary  times.  As  a  reminder,  L  is  the  beam  length  and 
v  is  the  particle  velocity  which  is  equal  to  PcQ.  The  tables  at  the  end  of  the  chapter  give 
the  times  in  normalized  times  and  are  indicated  by  Tj  (where  i  =  fa,b.c,d,l,2.3}). 

For  ease  of  reference,  the  above  limits  are  summarized  for  use  in  Appendix  A. 
Tables  1,  2,  and  3  from  (Ref.  2:  pp.  24,25,311,  give  the  expressions  for  z  \in  w  =  z  -  z) 
for  the  paths  to  the  right,  left,  and  centered  respectively.  The  expressions  for  a  and  b 
found  tn  Tabic  I.  are  given  :n  Equations  2.1b  and  2.1".  In  Taoie  2.  the  expressions  for 
aa  and  bb  are  given  in  Equations  2.18  and  2.19. 

The  magnetic  fields  can  now  be  calculated  for  any  location  P(z.s).  The  next 
chapter  uses  Appendix  A  to  generate  magnetic  field  (Cerenkov  pulse)  graphs  for  a 
couple  of  illustrative  examples.  Each  case  will  use  a  50  MeV  energy  level,  but  will  vary 
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Figure  2.11  Path  to  the  Right:  t  <  tfa. 

the  Au  L  values  in  order  to  study  any  differences.  Aiso,  the  developement  of  the  S-Z 
plane  is  easier  to  explain  since  establishing  the  origin  of  the  boundary  times  in 
Equations  2.12  through  2.15. 


Figure  2.16  Centered  About  the  Minimum. 
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Figure  2.P  tj  <  t  <  t,. 


Figure  2.18  t2  <  t  <  t,. 


TABLE  1 

LIMITS  OF  INTEGRATION:  PATH  TO  THE  RIGHT 


Tc  >  Tb: 
z'i 
zf 

Tc  <  Tb: 

z' 


‘f 

Tc  =  Tb: 
z . 


T  <  T  <  Th 

a  o 

0 


Ta  <  T  <  Tc 

0 


T  c  T  <  T. 

a  b 

0 


T.  <  T  <  T 

0  C 

b 

a 


TC  "  T  <  Tb 
0 


Tb  <  T  <  Td 

b 

L 


T  <  T  <  T.. 

C  U 

b 

L 

Tb  <  T  <  Td 

b 


TABLE  2 

LIMITS  OF  INTEGRATION:  PATH  TO  THE  LEFT 


T  <  T.: 

a  a 

z'i 

zr 

T.  <  T  <  T 

o  a 

aa 

L 

T  <  T  <  T. 

a  a 

0 

L 

T,  <  T  <  T. 

d  0 

0 

bb 

Ta  >  Td: 
z'. 

zf 

Tc  <  T  <  Td 
aa 

L 

Td  <  T  <  Ta 
aa 

bb 

Ta  <  T  <  Tb 
0 

bb 

Ta  =  Tb: 

z'i 

zf 

Tc  <  T  <  Ta 

aa 

L 

Td  <  T  <  Tb 

0 

bb 

TABLE  3 

LIMITS  OF  INTEGRATION:  CENTERED 

T!<T<T2 

z'.  ^  0 

z'_  >  0  z'+  ^  L  z'+  <  L 

z'i 

0 

z'_ 

z'f 

L  z'  + 

t2<t<t3 

First  Integral 

z'i 

0 

Z'.(Uj) 

zf  =  z'.(u2) 
Second  Integral 
z'i  =  z'_(u2) 
z'f 

L  z'  +  Iu^ 

III.  ANALYSIS  OF  THE  S-Z  PLANE 


A.  S/L  VERSES  Z/L 

After  establishing  the  pertinent  variables,  relative  observing  positions,  and  the 
different  situations  for  each  observing  position  in  the  previous  chapter,  the  S/L  vs.  Z/L 
graph  can  now  be  constructed  and  studied.  This  particular  graph  gives  some  insight  to 
the  integration  procedures  used  in  the  previous  chapter  since  the  limits  of  integration 
change  at  certain  times.  The  S/L  vs.  Z/L  graph  plots  out  the  conditions  where  the 
previously  derived  boundary  times  are  equal  and  also  presents  the  plane  of  radiation 
and  the  beam.  Figure  3.1  (page  36)  illustrates  this  graph  and  was  produced  using  the 
following  values:  J3'  =  1.1,  n  =  1.100057203,  L  =  150.0  meters,  and  Au/L  =  0.05. 
These  particular  values  spread  the  plotted  lines  out  and  allows  us  to  study  their  true 
shapes.  The  functions  that  plot  these  lines  are  obtained  by  setting  tb  =  t  (solid  line), 
ta  =  td  (dashed  line),  t,  =  ta  or  tb  =  td  (chain-dashed  line).  The  FORTRAN  program 
in  Appendix  B  creates  the  data  that  plots  out  each  function  in  Figure  3.1.  The 
conditions  that  must  exist  in  order  to  plot  the  above  functions  are  expressed  in 
Equations  3.1  to  3.3  for  tb  =  tc,  tfl  =  td,  and  ta  =  tc  respectively. 

Au  =  L  4-  P'[(s2  +  (z  -  L)2)'"'2  -  (s2  4-  z2)'/l]  (tb  =  tc)  (3.1) 

Au  -  P'[(s2  4-  z2)l/2  -  (s2  4-  (Z  -  L)2)’/2J  -  L  (ta  =  td)  (3.2) 

L  =  p'[(s2  +  z2)Vl  -  (s2  +  (z  -  L)V/21  (ta  =  tc)  (3.3) 

Equation  3.3  can  also  be  derived  by  setting  tb=  td  and  is  not  dependent  on  Au.  Figure 
3.2  :s  the  expanded  version  of  Figure  3.1  anu  is  shown  si  order  illustrate  where  remote 
ooserving  positions  fail. 

To  see  where  the  Cerenkov  region  (marked  center)  falls  on  the  S-Z  plane,  the 
value  of  0C  is  measured  from  the  horizontal  at  the  beam  discharge  point  (S/L  =  0,  Z/L 
=  0),  and  the  end  of  the  beam  (S/L  =  0,  Z/L  =  1).  By  projecting  a  line  from  each 
point  at  the  angle  0c  to  the  end  of  the  S-Z  plane  gives  the  Cerenkov  region.  In  Figures 
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3.1  and  3.2  this  region  falls  between  the  two  dotted  parallel  lines.  The  three  possible 
observer  regions  are  shown  in  Figure  3.1.  From  the  beam  discharge  point  to  the  first 
encountered  Cerenkov  boundary,  are  all  points  to  the  right  of  the  path.  In  the 
Cerenkov  region,  all  points  are  centered  about  the  minimum.  From  the  Cerenkov 
region  in  the  increasing  Z/L  direction  are  all  points  which  fall  to  the  left  of  the  path. 

B.  HYPERBOLAS  IN  THE  S-Z  PLANE 

Inspection  of  the  curves  in  Figures  3.1  and  3.2  leads  one  to  inquire  if  they  have  a 
simple  mathematical  expression.  As  mentioned  in  the  previous  section,  particular 
boundary  times  were  set  equal  in  order  to  obtain  these  three  functions  of  interest.  The 
initial  mathmatical  expression  that  appears  after  equating  these  times  is  given  by 
Equation  3.4,  where  y  =  S/L  and  x  =  Z/L. 

(x2  +  y2)'/2  -  ((x  -  l)2  +  y2)‘/2  =  Q  (3.4) 

The  value  of  Q  is  unique  to  each  function  and  is  always  a  positive  constant.  Equations 
3.5,  3.6,  and  3.7  are  the  three  expressions  for  Q  that  originate  when  setting  tfe  =  t 
(right),  ta  =  t .  (left),  and  ta  =  t  (center)  respectively. 

Qr  -  [1  -  Au/L]/p'  (3.5) 


Q1  =  [1  4-  Au/L]/P' 


Qc  =  1/P  (3.7) 

Equation  3.4  is  the  difference  between  two  distances  equated  to  a  constant.  This 
constitutes  the  standard  form  of  a  hyperbola.  In  order  to  write  a  program  for  plotting 

hyperbolas,  an  exDlicit  expression  for  y  as  a  function  of  x  is  needed.  After  some 
cumbersome  algebra.  Equation  3.4  becomes  Equation  3.3. 

y  =  l((x2  -  (x-  l)2  -  Q2)2/4Q2)  -  (x  -  l)2]"2  (3.8) 
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The  next  step  is  to  take  Equation  3.8  and  work  it  into  the  classical  form  of  a  hyperbola 
where  there  is  an  axis  translation.  This  allows  one  to  find  the  coordinates  (x0,y  )  of  the 
axis  translation.  Equation  3.9  gives  the  classical  form  where  the  hyperbola  is 
symmetric  about  xQ  and  y  . 

(x-x0)2/a2  -  (y  -  y0)2/b2  =  i  (3.9) 

The  best  approach  to  achieving  the  classical  form  is  to  work  with  Equation  3.8. 
By  carrying  out  the  square  in  the  (x  -  l)2  term  and  performing  other  algebraic 
operations,  a  general  form  for  conical  sections  is  recognized.  Equation  3.10  takes  on 
the  form  Ax2  +  Bx  +  Cy2  +  Dy  +  E  =  0. 


yz  -  (1/QZ  -  l)xz  +  [2  -  (1  +  Qz)/Qz]x  -  (1  -  2 Qz  +  Q4)/ 4Q2  =  0 


(3.10) 


A  hyperbola  exists  if  AC  is  negative.  An  ellipse  or  parabola  exists  if  AC  is  positive  or 
zero  respectively.  From  Equation  3.4,  the  only  possible  shapes  that  form  on  the  S-Z 
plane  are  hyperbolas.  To  determine  the  conditions  for  the  existence  of  a  hyperbola,  a 
closer  look  at  the  coefficient  of  the  x2  term  is  required. 

Since  the  only  true  shapes  on  the  S-Z  plane  are  hyperbolas,  then  the  coefficient 
of  x2  must  be  positive.  In  other  words  the  inequality  3.11  must  be  satisfied. 


1/Q  >  1 


(3.11) 


Using  the  expressions  for  Q  in  Equations  3.5,  3.6,  and  3.7  with  the  particular  values  of 
Au,  L,  and  P',  the  inequality  is  easier  to  satisfied  when  considering  Equations  3.5  and 
3.7.  However,  when  using  Equation  3.6,  difficulties  are  more  likely  to  arise  in 
obtaining  the  expected  hyperbolas.  Substituting  in  Q1  for  Q  in  the  inequality  3.11, 
allows  one  to  see  how  the  values  of  P',  Au,  and  L  govern  the  inequality.  The  inequality 
3.12  is  'he  result  of  this  substitution. 


P'  >  1  +  Au/L 


(3.12) 


If  the  inequality  3.12  is  satisfied,  the  time  equalities  which  give  rise  to  Equations  3.1  to 
3.3  are  satisfied.  The  most  common  values  encountered  with  the  variables  P'  and  Au  L 
cause  the  inequality  3.12  to  be  unsatisfied. 
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In  summary,  hyperbolic  shapes  exist  in  the  S-Z  plane  when  the  above  inequalities 
are  satisfied.  Without  explicitly  defining  the  values  of  a2  and  o2,  it  can  easily  be  shown 
from  Equation  3.10  that  the  hyperbolas  take  on  the  form 

tx  -  0.5)2/a2  -  v2b2  =  1.  t3.15) 

The  values  of  a  and  b  depend  on  the  function  that  is  being  plotted  and  are 
subsequently  functions  of  the  corresponding  value  of  Q. 

C.  CERENKOV  PULSE  SHAPES 

From  Appendix  A,  Cerenkov  pulses  are  generated  from  knowing  many  variables, 
particularly  the  location  of  the  observer  with  respect  to  the  beam  discharge  point. 
From  here,  the  relationship  between  the  Cerenkov  pulse  shape  and  position  relative  to 
time  boundaries  in  the  S-Z  plane  can  be  studied  graphically.  The  best  way  to  study  this 
relationship  is  to  methodically  choose  points  on  the  S-Z  plane  that  cross  the  critical 
time  and  position  boundaries.  Other  points  are  also  chosen  in  order  to  aid  in 
developing  general  trends  in  the  way  the  pulse  shape  changes  as  the  observing  points 
move  across  the  S-Z  plane.  The  plus  signs  in  Figures  3.1  and  3.2  indicate  the  points 
where  the  Cerenkov  pulses  are  calculated.  Each  observing  position  is  numbered  in 
order  to  simplify  the  discussion  in  the  conclusions.  As  indicated  in  the  right  margin,  p' 
-  1.1,  Au.'L  =  DU/L  =  0.05,  and  L  =  150.0  meters.  Choosing  a  constant  S/L  value 
of  0.26  and  different  Z/L  values  near  a  critical  boundary,  gives  a  good  representation  of 
how  the  pulse  shape  varies.  Figures  3.3,  3.4,  and  3.5  fall  behind  the  beam  discharge 
point.  The  magnitudes  have  been  scaled  up  by  a  factor  of  ten  in  order  to  better 
examine  the  shape  of  the  pulse.  Figures  3.6  to  3.13  continue  to  move  in  the  plus  Z/L 
direction  while  crossing  critical  boundaries.  Since  the  magnitude  of  the  Cerenkov  pulses 
are  large  enough,  their  values  are  not  scaled  up.  By  visual  inspection,  one  can  observe 
that  the  pulse  peaks  increase  as  you  approach  the  Cerenkov  region.  After  the 
maximum  peak  occurs,  the  magnitude  of  the  peaks  decrease  with  increasing  Z/L.  The 
jr served  oreuxs  in  the  pulse  shapes  occur  at  the  deiined  normalized  boundary  times 
T_,  T. ,  T,,  ana  Tr 

aT  b  c’  a 

The  remote  observing  points  are  shown  in  Figures  3.14  to  3.20.  These  points  are 
indicated  on  Figure  3.2  and  are  selected  to  fall  around  critical  boundaries.  Figures  3.14 
and  3.15  are  located  to  the  right  of  the  path.  The  pulse  shapes  peak  at  Tb  then 


gradually  decreases  in  magnitude  to  Tc  before  decreasing  to  zero  at  Td.  Figures  3.18 
and  3.19  are  points  which  fall  in  the  Cerenkov  region.  The  pulses  begin  at  time  Tj  and 
end  at  time  T3  with  a  peak  at  T2.  The  peaks  of  each  of  the  pulses  appear  flat,  and 
show  no  gradual  decline  as  did  the  pulses  to  the  right  of  the  path.  Figure  3.20  is  a 

pulse  ;c  the  left  of  the  path  with  shape  characteristics  that  are  similar  to  Figures  3. IS 

and  3.19. 

To  continue  this  study  of  the  relationship  between  the  S-Z  plane  and  the 
Cerenkov  pulse  shapes,  the  parameter  Au/L  (DU/L)  is  changed  from  0.05  to  0.15. 
Figures  3.21  and  3.22  are  shown  as  were  Figures  3.1  and  3.2  with  the  regions  marked 
and  the  points  (plus  signs)  for  evaluation  indicated.  However,  the  dashed  line  (where  t 
=  td)  is  not  present.  This  is  due  to  the  fact  that  the  inequality  3.12  is  not  satisfied  as  a 
result  of  the  new  Au/L  ratio.  Thus,  only  two  hyperbolic  shapes  exist  in  Figures  3.21 
and  3.22. 

Figures  3.23  to  3.34  are  the  same  in  shape  as  Figures  3.3  and  3.13.  The 
magnitudes  are  larger  even  though  the  points  from  which  these  pulses  were  generated 
are  a  longer  distance  away  from  the  discharge  point.  Figures  3.35  through  3.41  are  the 
pulses  at  the  points  indicated  on  Figure  3.22.  The  shapes  are  again  somewhat  expected 
for  each  case.  The  points  located  to  the  right  of  the  path  peak  and  break  at  the  critical 
times,  as  does  the  points  located  to  the  left  and  on  the  center.  However,  the  pulses 
located  in  the  Cerenkov  region  and  to  the  left  of  the  path  have  flattop  peaks.  As  an 
obvious  observation,  the  remote  points  in  each  Au/L  case  have  later  times  at  which  the 
pulse  is  starting  to  be  observed.  This  is  due  to  the  longer  distances  that  the  emissions 
have  to  travel  in  order  to  get  to  the  observing  positions. 

The  next  approach  which  aids  in  gaining  better  insight  in  the  study  of  the  pulse 
shapes,  is  to  examine  points  that  are  an  equal  distance  from  the  origin  of  the  beam. 
Figures  3.42  and  3.43  originate  from  generating  the  S-Z  plane  using  the  Au/L  ratio  of 
0.05  as  in  Figures  3.1  and  3.2.  The  reason  for  using  this  case  as  an  example  is  because 
all  three  hyperbolic  shapes  exist.  Figures  3.42  and  3.43  represent  the  close  and  far 
observing  positions  -espectivelv.  The  exoressions  close  and  far  observing  positions  are 


:n  terms  af  their  relative  magnitudes.  Close  simolv  describes 


ooservini 


position  that  is  nearer  to  the  beam.  The  points  where  the  pulses  are  calculated  are 
again  indicated  by  a  plus  sign.  The  quarter  circle  (solid  line)  has  a  radius  of  curvature 
pointed  toward  the  beam  discharge  point.  Figures  3.44  to  3.48  are  the  pulse  shapes 
from  Figure  3.42.  Figures  3.44  and  3.45  fall  to  the  right  of  the  path.  Figures  3.46  and 


3.47  lie  in  the  Cerenkov  region  (centered),  and  Figure  3.4S  falls  to  the  left  of  the  path. 
Figures  3.49  to  3.57  correspond  to  the  points  on  Figure  3.43.  The  arrival  time  of  the 
emitted  signal  is  the  start  of  each  pulse  and  is  the  same  when  comparing  two  points  in 
the  same  region.  The  earliest  arrival  times  occur  to  the  left  of  the  path,  when 
considering  the  equal  distance  criterion,  while  the  latest  arrival  times  occur  to  the  right 
of  the  path.  As  a  quick  observation,  the  magnitudes  of  the  pulses  get  larger  and  the 
pulse  widths  narrower  as  the  Cerenkov  region  is  approached.  More  detailed 
conclusions  and  observations  will  follow  in  the  next  chapter. 
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Figure  3.8  From  Figure  3.1:  Pulse  6. 
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Figure  3.9  From  Figure  3.1:  Pulse  7. 
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Figure  3.10  From  Figure  3.1:  Pulse  8. 
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Figure  3.14  From  Figure  3.2: 
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Figure  3.16  From  Figure  3.2:  Pulse  14. 
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Fisiure  3.17  From  Fisure  3.2:  Pulse  15. 
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Figure  3.20  From  Figure  3.2:  Pulse  18 
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Figure  3.25  From  Figure  3.21:  Pulse  21. 
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Figure  3.26  From  Figure  3.21:  Pulse  22. 
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Figure  3.28  From  Figure  3.21:  Pulse  24. 
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Figure  3.30  From  Figure  3.21:  Pulse  26. 


5 


1.0  1.5  2.0 

TIME  (NORMALIZED) 


Ficure  3.31  From  Ficurc  3.21: 
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Figure  3.46  From  Figure  3.42:  Pulse  40. 
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Figure  3.49  From  Figure  3.43:  Pulse  43. 
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Figure  3.54  From  Figure  3.43:  Pulse  48. 
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Figure  3.56  From  Figure  3.43:  Pulse  50. 
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v.  3aiaOROIND 

The  shape  of  the  B  held  t  Cerenkov  pui>c-  ,>  known  to  he  depeude:::  or.  On 

position  relative  to  the  mimmum  of  the  function  utz  >  However,  based  on  the 

conclusions  in  IRef.  2  p  44],  it  was  suspected  ’.hat  the  pulse  shape  ;s  also  related  ’o  its 

position  relative  to  the  time  boundary  regions.  Again,  these  regions  arc  where  t.  =  t  . 

t  *  t,,  and  t  =*  t  The  recommendations  that  followed  from  this  suspicion  were  to 
a  d  c  a 

use  different  values  for  Au  L  and  to  generate  Cerenkov  pulses  at  points  close  to  these 
time  boundaries.  In  the  previous  chapter.  Figures  3.1  to  3.57  allows  one  to  compare 
various  shapes  for  different  situations  and  Au  L  values.  These  figures  also  aid  m 
drawing  conclusions  about  the  pulse  shapes  and  sizes. 

B.  OBSERVATIONS 

1.  General  Trends 

So  far  only  a  few  brief  remarks  have  been  made  with  regard  to  the  shapes  and 
magnitudes  of  the  Cerenkov  pulses  in  the  previous  chapter.  Here,  more  specific 
conclusions  are  drawn  and  some  general  observations  are  discussed.  To  make  this 
easier,  a  useful  way  to  study  and  compare  the  pulse  shapes,  and  at  the  same  time  see 
where  they  fall  with  respect  to  the  critical  boundaries  is  to  use  Figures  4.1  to  4.4.  The 
pulses  are  numbered  in  ascending  order  from  left  to  right  starting  with  the  position 
closest  to  the  origin  and  then  moving  out  to  greater  distances.  Figures  4.1  and  4.2 
represent  the  close  and  far  field  regions  respectively  and  correspond  to  the  Au/L  value 
of  0.05.  They  are  simply  Figures  3.1  and  3.2  with  the  Cerenkov  pulses  superimposed. 
The  pulses  are  sketched  to  show  relative  magnitudes,  shapes,  and  pulse  widths.  Figures 
4.3  and  4.4  correspond  to  Figures  3.21  and  3.22.  The  pulse  numbering  continues  where 
the  last  pulse  number  stops  in  Figure  4,2. 

Starting  with  the  puise  on  the  far  left  (number  1)  in  Figure  4.1  and  proceeding 
in  the  positive  Z/L  direction,  the  pulses  have  the  same  basic  shape.  This  means  that  the 
magnitude  of  the  pulse  goes  from  zero  to  a  peak  value,  then  declines  to  a  smaller  value 
before  falling  off  to  zero.  Pulses  1  to  5  fall  to  the  right  of  the  path  and  do  not  cross  a 
critical  boundary.  Each  pulse  begins  at  time  ta  and  then  increases  to  a  peak  value  at  tb> 
From  here,  the  pulse  declines  in  a  non-linear  fashion  to  a  smaller  value  at  t  before 
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reiav.se  increase  an-  starts  to  Mend  into  the  decline  from  the  peak,  magnitude  to  the 
end  of  the  puNe  Since  the  observing  position  is  approaching  the  tb  *  t  boundary 
■s .me  .  me  .:. teg:  it.cn  pro-ess  ta^cs  t..e  . a.ue  o;  tne  integral  to  uie  peux  at  t.  and 
tneti  ba-k  down  to  zero  For  referen-e.  puise  "  corresponds  to  Figure  3.4.  Crossing  the 
*b  ~  C  boundary  to  pulse  P,  shows  that  the  t.  value  has  vanished  and  the  shape  takes 
the  form  of  a  spike.  Pulses  4.  10,  and  i  1  on  Figure  4.1,  cross  two  time  boundaries 
before  leaving  the  Cerenkov  region.  The  last  pulse  falls  to  the  left  of  the  path.  There  is 
no  cnange  in  snape;  the  last  three  pulses  are  spikes.  Their  magnitudes  decrease  and  the 
pulse  widths  increase  with  increasing  Z  L.  Pulses  S  to  1 1  correspond  to  Figures  3.10  to 
31' 


The  above  discussion  considered  Figure  4. 1  and  moved  from  left  to  right 
examining  each  pulse.  Comments  were  made  when  a  critical  boundary  was  crossed  and 
two  pulses  were  compared  when  they  were  positioned  on  opposite  sides  of  a  boundary. 
This  procedure  for  examining  Figure  4.1  w’as  necessary  but,  will  not  be  repeated  for 
Figures  4.2  to  4.4.  Instead,  only  the  general  trends  in  shapes,  magnitudes,  and  pulse 
widths  will  be  noted. 

Figure  4.2  illustrates  the  pulses  that  were  generated  at  larger  distances  from 
the  beam  discharge  point.  The  maximum  value  of  the  pulse  increases  from  left  to  right 
up  to  pulse  number  15  (Figure  3.17),  and  from  there  a  decline  begins.  The  pulse 
widths  decrease  and  then  increase  in  the  same  manner  that  the  magnitudes  increase  and 
decrease.  The  shapes  are  similar  to  the  close  field  observations  with  the  exception  of 
one  minor  difference.  The  smaller  value  at  tc  is  almost  as  large  as  the  peak  magnitude 
of  the  pulse.  After  ’he  slightly  rounded  spike  puise  <  !5i  occurs,  the  remaining  pulses 
ore  flattops  is  one  proceeds  to  larger  distances.  From  the  maximum  peak  magnitude, 
the  rest  of  the  shapes  show  decreasing  magnitudes  and  increasing  pulse  widths  with 
increasing  Z/L. 

Figures  4.3  and  4.4  also  show  the  close  and  far  fields  respectively.  They  are 
similar  to  Figures  4.1  and  4.2  with  the  exception  that  the  value  for  Au/L  is  0.15  vice 
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0.05.  Again,  the  general  trends  in  magnitudes  and  pulse  widths  are  the  same.  However, 
the  pulses  further  to  the  right  take  on  flattop  peaks  where  spikes  were  observed  in  the 
previous  two  figures.  Other  than  this,  analyzing  Figures  4.3  and  4.4  would  be  repeating 
what  was  done  earlier  in  Figures  4.1  and  4.2. 

Z.  Specific  Analysis 

At  this  point,  observations  have  been  made  by  holding  S/L  constant  and 
varying  Z/L.  This  means  that  the  distance  from  the  beam  discharge  point  is  changing 
from  one  pulse  to  the  next.  Figures  3.42  and  3.43  illustrate  the  close  and  far  observing 
positions  respectively  where  the  distance  from  the  discharge  point  is  constant.  Inis 
allows  one  to  see  how  the  pulses  change  with  respect  to  location  in  the  S-Z  plane 
without  the  effects  of  a  varying  range  from  the  emission  source.  For  reference,  Figures 
3.44  to  3.48  correspond  to  pulses  38  to  42  in  the  close  positions,  while  Figures  3.49  to 
3.57  correspond  to  pulses  43  to  51  in  the  far  positions. 

The  pulses  from  Figures  3.42  and  3.43  yield  no  new  information  with  regard  to 
shape,  but  it  is  obvious  that  the  peak  magnitude  occurs  in  the  middle  of  the  Cerenkov 
region  as  does  the  minimum  pulse  width.  Without  the  aid  of  the  equal  distance 
analysis,  this  observation  would  be  only  a  safe  assumption.  Figures  4.5  and  4.6  are 
plots  of  the  peak  magnitude  and  pulse  width  verses  the  angle  measured  from  the 
vertical  axis  S/L.  They  both  show  that  the  angle  at  which  the  peak  magnitude  and 
minimum  pulse  width  occur  is  the  compliment  of  the  Cerenkov  angle.  In  other  words, 
measuring  the  angle  from  the  horizontal  axis,  would  yield  the  peak  magnitude  and  the 
minimum  pulse  width  at  the  Cerenkov  angle.  This  is  consistent  with  Cerenkov 
radiation  theory. 

C.  SUMMARY 

One  conclusion  in  [Ref.  2:  p.  45],  that  can  be  reinforced  after  analyzing  a.  ■ 
pulse  shapes  in  this  paper,  is  that  a  smaller  Au.'L  value  will  yield  more  m — .  •  - 
shapes.  Still,  there  are  some  other  points  that  need  to  be  brought  out  m  - 

pointed  pulse  shapes,  fall  closer  to  a  time  bcundarv  ■han  -he  ~  .  - 

S-Z  plane.  Figures  3. in  puisc  ■  *.  ' . 5* >  -ru.r-e  >  .  '  -  .  .  - 

are  specific  examples  of  tins  observation  a  no  appi.es  :  ,o. 
exception  to  some  of  the  cases  involving  the  t  -  t  w  >.•  : 
sense  because  on  the  boundary  there  are  .  >  :w  a-  .n- 
through  during  the  integration  prtxc".  1  -  ew 
rise  to  a  peak  at  t  from,  t  and  ♦hen  d-„\  :  ••  • 

A  o 
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move  away  from  the  time  boundary,  there  is  enough  distinction  between  the  two  times 
for  the  computer  to  notice  two  distinct  magnitudes.  This  explains  why  the  value  at  tc> 
continues  to  approach  the  magnitude  of  tfe,  when  the  observing  positions  change  from 
left  to  right  on  the  S-Z  plane  when  starting  to  the  right  of  the  path.  In  other  words,  if 
an  observing  position  falls  on  a  time  boundary  where  two  critical  times  are  equal,  then 
there  exists  a  unique  value  for  the  magnitude  of  the  pulse.  In  short,  all  this  means  that 
the  general  shape  of  the  Cerenkov  pulse  depends  on  the  location  of  the  observer 
relative  to  a  time  boundary.  In  particular,  for  positions  near  or  on  a  time  boundary, 
the  observed  pulses  are  more  spike-like.  However,  the  magnitude  and  the  pulse  width 
are  dependent  on  the  observer's  position  relative  to  the  Cerenkov  region. 
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APPENDIX  A 

CERENKOV  PULSE  PROGRAM 


*THIS  PROGRAM  IS  LOCATED  ON  THE  NPS  MAIN  FRAME  UNDER  THE  FILENAME  * 
* 1 CERENKV '  AND  IS  WRITTEN  USING  WATFOR  77.  AFTER  ENTERING  CERTAIN  * 
*VALUES ,  DATA  FOR  THE  MAGNITUDE  OF  THE  CERENKOV  PULSE  IS  GENERATED* 
*AND  IS  THEN  PLOTTED  BY  ANY  PLOTTING  ROUTINE.  THESE  MAGNITUDES  ARE* 
*PLOTTED  AGAINST  NORMALIZED  TIME  AND  WHEN  MULTIPLIED  BY  A  CONSTANT* 
*YOU  CAN  GET  ACTUAL  B-FIELD  OR  E-FIELD.  * 

DOUBLE  PRECISION  N,  U1 ,  U2,  BETA,  CO,  ROE,  A,  G,  CE,  BPRME ,  POSZ 
DOUBLE  PRECISION  Rl,  R2 ,  Al,  RSETME ,  SRTE1 ,  SRTE2 ,  A2,  D,  E,  DD 
DOUBLE  PRECISION  EE,  Q,  TA,  TB,  TC,  TD,  F,  ZPI ,  ZPF,  WPI ,  WPF 
DOUBLE  PRECISION  W1 ,  W2 ,  ZC,  RC,  ZPC,  T1 ,  T2,  T3,  DELR,  ZPM,  B1 
DOUBLE  PRECISION  B2,  El,  E2 ,  SI,  XX,  YY,  S,  L,  Z,  TPRME ,  B,  BMAX 
DOUBLE  PRECISION  TMAX,  TMIN ,  GPRME ,  THETA1 ,  THETA2 ,  TINC,  TAU 
DIMENSION  TPRME (9000),  B(9000) 

INTEGER  I,  J,  IMAX 

*****A*********K?t>i**j:******************************************x**jt 

*  ALL  THE  PERTINENT  VARIABLES  ARE  ENTERED  HERE  AND  SELF-EXPLAINS  * 

*  THE  NEED  FOR  THE  VARIABLE  WITH  TERMINAL  INTERACTION.  ONE  POINT  * 

*  THAT  NEEDS  TO  BE  CLARIFIED,  IS  HOW  TO  ENTER  U1  AND  U2.  U1  WILL  * 

*  BE  ZERO  AND  U2  WILL  BE  A  NEGATIVE  VALUE.  * 

PRINT  *,  'ENTER  N1 
READ  *  N 

PRINT  *,  ‘ENTER  ROE* 

READ  *  ROE 

PRINT  *,  'ENTER  BETA' 

READ  *  BETA 

PRINT  *,  'ENTER  BEAMLENGTH' 

READ  *  L 

PRINT  *,  'ENTER  U1  SHOULD  USUALLY  BE  ZERO' 

READ  *,  U1 

PRINT  *,  'ENTER  U2' 

READ  *,  U2 

PRINT  *,  'ENTER  S  POSITION' 

READ  *  S 

PRINT  *,  'ENTER  Z  POSITION' 

READ  *  Z 

PRINT  *,  'ENTER  CO  IN  M/NSEC' 

READ  *,  CO 

V  =  BETA*CO 

BPRME  =  N*BETA 

Q  =  ( BPRME **2 .0)  -  1.0 

Rl  =  SQRT( (S**2 .0)  +  (Z**2 . 0) } 

R2  =  SQRT((S**2.0)  +  (Z  -  L)**2.0) 

si  =  (§**2.o)*q 

TAU  =  L/7 

*THE  30UNDARY  TIMES  ARE  COMPUTED  AND  USED  THROUGHOUT  THE  PROGRAM  * 
*DURING  THE  INTEGRATION  PROCESS  FOR  CALCULATING  THE  MAGNITUDE  OF  * 
*THE  CERENKOV  PULSE.  THE  WRITE  STATEMENT  DISPLAYS  THE  NORMALIZED  * 
*TIMES  WHICH  SERVES  AS  AN  AID  WHEN  USING  VARIOUS  GRAPHICS  ROUTINES* 
*FOR  PLOTTING  THE  PULSES.  * 

**********************************  *******************  A***********  *  * 

TA  =  ( (BPRME*R1  -  Ul)/V) 

TB  =  ( (BPRME*R1  -  U2)/V) 

TC  =  ((L  +  (BPRME*R2 )  -  Ul)/V) 

TD  =  ((L  +  (BPRME*R2)  -  U2)/V) 


WRITE (6 ,110)  TA/TAU,  TB/TAU,  TC/TAU,  TD/TAU 

110  FORMAT(4F10.5) 

£  **** A******* A*********************************************** ******* 

C  INITIALIZE  THE  DIMENSIONED  VARIABLES  FOR  TIME  AND  MAGNITUDE  -  B  * 

Q  ************************************************************ ******* 

DO  10  I  =  1,200 
TPRME(I)  =  0.0 
3(1)  =  0.0 
10  CONTINUE 

C  ******************************************************************* 

C  *THE  CERENKOV  ANGLE  CE  IS  CALCULATED  AND  BASED  ON  THE  OBSERVING  * 
C  *POSITIONS,  THE  VALUES  OF  THETA1  AND  THETA2  ARE  CALCULATED.  FROM  * 
C  *HERE  WE  CHECK  TO  SEE  IF  WE  ARE  TO  THE  LEFT,  RIGHT,  OR  CENTER.  * 

Q  ******************************************************************* 

CE  =  DACOS ( 1 . 0/BPRME ) 

IF  (Z  .LE.  0.0)  THEN 
POSZ  =  ABS(Z) 

A  =  POSZ/S 
THETA1  =  DAT AN (A) 

THETA1  =  THETA1  +  1.570796327 
ELSE 

A  =  S/Z 

THETA1  =  DAT AN (A) 

END  IF 

IF  (Z  .LT.  L)  THEN 
GPRME  =  (L  -  Z)/S 

THETA2  =  DATAN ( GPRME )  +  1.570796327 
END  IF 

IF  (Z  .EQ.  L)  THEN 

THETA2  =  1.570796327 
END  IF 

IF  (Z  .GT.  L)  THEN 
GPRME  =  S/(Z  -  L) 

THETA2  =  DATAN (GPRME) 

END  IF 

C  ******************************************************************* 

C  *THE  CERENKOV  ANGLE  IS  COMPARED  TO  THETA1 .  THE  IF  STATEMENT  CHECKS* 
C  *TO  SEE  IF  WE  ARE  TO  THE  RIGHT  AND  THUS  DISPLAYS  ON  THE  SCREEN.  * 

C  ******************************************************************* 

IF  (THETA1  .GT.  CE)  THEN 
WRITE (6 , 111 ) 

111  FORMAT  ('PATH  TO  THE  RIGHT') 

TINC  =  (TD  -  TA)/200 .0 
DO  20  I  =  1,200 

IF  (I  .EQ.  1.0)  THEN 
TPRME(I)  =  TA 
B ( I )  =  0.0 
ELSE 

TPRME(I)  =  TA  +  (REAL(I)*TINC) 

END  IF 

IF  (TPRME(I)  .GE.  TD)  THEN 
B(I)  =  0.0 
GO  TO  100 
END  IF 

Ai  =  TJ1  *  7*T?RME  ( I ) 

A 2  =  'J2  +  7*T?RME  ( I ; 

D  =  ( (BPRME*X2 . 0 )*Z)  -  Al 
DD  =  ( (BPRME**2 .0)*Z)  -  A2 
El  =  ((Z  -  Al)**2.0)  -  SI 
E2  =  ((Z  -  A2 ) **2 . 0 )  -  SI 
IF  (El  .LT.  0.0)  THEN 
El  =  0.0 
END  IF 

IF  (E2  .LT.  0.0)  THEN 
E2  =  0.0 
END  IF 

SRTE1  =  SQRT(El) 


E  =  BPRME^SRTEl 
EE  =  BPRME*SRTE2 
ZPI  =  (DD  +  EE)/Q 
ZPF  =  (D  +  E)/Q 
IF  (ZPI  .LT.  0.0  )  THEN 
ZPI  =0.0 
END  IF 

I?  (ZPF  .LZ.  0.0)  THEN 
ZPF  =  0.0 
END  IF 

***************************************** ************ ************** 
*WHEN  TO  THE  RIGHT  OF  THE  PATH  TC  >  TB  IS  ONE  OF  THE  CASES  WE  HAVE* 
*TO  CHECK  FOR  AND  THEN  CARRY  OUT  THE  INTEGRATION.  * 

****************************** ************** *********************** 


TA  .AND.  TPRME(I)  .LT.  TB)  THEN 


0) )*(YY-XX) 
TPRNE(I) 


LT.  TC)  THEN 


IF  {TC  .GT.  TB)  THEN 
IF  (TPRME(I)  .GT. 

WPI  =  Z  -  0.0 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN(W2 
B(I )  =  (R0E*N* (BETA**2 
END  IF 

IF  (TPRME(I)  .GT.  TB  .AND. 

WPI  =  Z  -  ZPI 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN(W2) 

3(1)  =  RCE’LN* ( SETA**2 . 0 )  * ( YY  -XX) 

END  IF 

IF  (TPRME(I)  .GT.  TC  .AND.  TPRME(I)  .LT.  TD)  THEN 
WPI  =  Z  -  ZPI 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN(W2) 

=  ROE*N*(BETA**2.0)*(YY  - 


END 
END  IF 


B(I) 

IF 


XX) 


******************************************************************* 

* ANOTHER  CASE  FOR  THE  PATH  TO  THE  RIGHT  IS  TC  LESS  THAN  TB.  * 

******************************************************************* 


IF  (TC  .LT.  TB)  THEN 

IF  (TPRME(I)  .GT.  TA  .AND.  TPRME(I)  .LT.  TC)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  OAT  AN  ■!  W2  ! 

B.'T)  =  SCZ^N* ' BETA^^Z . 0 ) * ' YY  -  XX) 

END  IF 

IF  (TPRME(I)  .GT.  TC  .AND.  TPRME(I)  -LT.  T3)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN(W2) 

B(I)  =  ROE*N*(BETA**2.0)*(YY  -  XX) 

END  IF 

IF  (TPRME(I)  .GT.  TB  .AND.  TPRME(I)  .LT.  TD)  THEN 
WPI  =  Z  -  ZPI 
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END  IF 
END  IF 


W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN{W2) 

B(I)  =  R0E*N*(BETA**2.0)*(YY 


^FINALLY  THE  THIRD  CASE  FOR  THE  PATH  TO  THE  RIGHT.  TC  =  TB  * 

'k-k*-k-k'k'k**-kj<ic'k-k*:-k-k1ck-k'k7**  ***  ******  *****************;*;*******  ****** 

IF  (TB  .EO.  TC)  THEN 

IF  (TPRME(I)  .GT.  TA  .AND.  TPRME(I)  .LT.  TB)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATANt W2) 

B(I)  =  ROE*N* ( BETA**2 . 0 ) * ( YY  -  XX) 

END  IF 

IF  (TPRME(I)  .GT.  TB  .AND.  TPRME(I)  .LT.  TD)  THEN 
WPI  =  Z  -  ZPI 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DA TAN ( W2 ) 

B(I)  =  R0E*N*(BETA**2.0)*(YY  -  XX) 

END  IF 
END  IF 

0  CONTINUE 

END  IF 

***************** A********* A* A* ********************* A****** ******** 

*THE  CERENKOV  ANGLE  IS  COMPARED  TO  THETA2  TO  SEE  IF  WE  ARE  TO  THE  * 
*LEFT  OF  THE  PATH.  THE  'WRITE'  WILL  INDICATE  THIS  ON  THE  SCREEN.  * 
******************************************************************* 

IF  (THETA2  .LT.  CE)  THEN 
WRITE (6 , 112) 

2  FORMAT ('PATH  TO  THE  LEFT') 

TINC  =  (TB  -  TC)/2Q0.0 

DO  25  I  =  1,200 

IF  (I  .EQ.  1.0)  THEN 
TPRME(I)  =  TC 
B ( I )  =  0.0 
ELSE 

TPRME(I)  =  TC  +  (REAL (I ^T INC) 

END  IF 

IF  (TPRME(I)  .GE.  TB)  THEN 
B(I)  =  0.0 
GO  TO  100 
END  IF 

A1  =  U1  +  (V*TPRME ( I ) ) 

a 2  =  j2  -  ,v'-t?rme(: ,  ) 

o  =  3?r:ie^-2  . ;  i  rz  •  -  ai 

ID  -  SPF.ME'^D  .  I  ;  )  -  A2 

II  =  . ,Z  -  k2)r~Z.j.  ~  S. 

E2  =  ((Z  -  A2 ) **2 . 0 )  -  SI 
IF  (El  .LT.  0.0)  THEN 
El  =  0.0 
END  IF 

IF  (E2  .LT.  0.0)  THEN 
E2  =  0.0 
END  IF 

SRTE1  =  SQRT(El) 

SRTE2  =  SORT (E2) 

E  =  BPRME*SRTE1 
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EE  =  BPRME*SRTE2 
ZPI  =  (D  -  E)/Q 
ZPF  =  (DD  -  EEJ/Q 
IF  (ZPI  .LT.  0.0)  THEN 
ZPI  =  0.0 
END  IF 

IF  (ZPF  .LT.  0.0)  THEN 
ZPF  =  0.0 
END  IF 

****************************  A**********  *******************  *  ******** 

* CHECK  FOR  ONE  OF  THREE  CASES  THAT  ONE  CAN  RUNINTO  WHEN  TO  THE  * 
*LEFT  OF  THE  PATH.  THIS  PARTICULAR  CASE  IS  TA  <  TD.  * 

********************************************  *************  ****  **  **** 

IF  (TA  .LT.  TD)  THEN 

IF  ( TPRME ( I )  .GT.  TC  .AND.  TPRME(I)  .LT.  TA)  THEN 
WPI  =  Z  -  ZPI 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN (W2) 

B(I)  =  ROE^N*(BETA**2.0)*(YY  -  XX) 

END  IF 

IF  (TPRME(I)  .GT.  TA  .AND.  TPRME ( I )  .LT.  TD)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN ( W2 ) 

3(1)  =  RCEKN* (BETA**2 . 0 ) * ( YY  -  XX) 


END  IF 

IF  (TPRME ( I )  .GT.  TD  .AND.  TPRME ( I )  .LT.  TB)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN (Wl) 

XX  =  DATAN  W2) 

B(I)  =  ROE*N*(BETA**2.0)*(YY  -  XX) 

END  IF 
END  IF 

******************************************************************* 

*THE  SECOND  CASE  FOR  THE  PATH  TO  THE  LEFT  IS  TA  >  TD.  * 

******************************************************************* 

IF  (TA  .GT.  TD)  THEN 

IF  (TPRME ( I )  .GT.  TC  .AND.  TPRME ( I )  .LT.  TD)  THEN 
WPI  =  Z  -  ZPI 
WPF  =  Z  -  L 
Wl  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN (Wl) 

XX  =  DATAN  W2) 

3a)  =  ROEai*(BETA**2.0)*(YY  -  XX) 

END  IF 

IF  TPRME  I'.  .FT.  TD  .AND.  TPRME  a '  .IT.  TA  i  THEN 
WPI  =  Z  -  DPI 
WPF  =  Z  -  ZPF 
Wl  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN (W2) 

B a )  =  ROE*N*(BETA**2 . 0)*(YY  -  XX) 

END  IF 

IF  (TPRME ( I )  .GT.  TA  .AND.  TPRME ( I )  .LT.  TB)  THEN 


W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN ( W2 ) 

B(I)  =  ROE*N*(BETA**2.0)*(YY  -  XX) 

END  IF 
END  IF 

******************************************* *********** ******* ****** 
*THE  THIRD  CASE  FOR  THE  PATH  TO  THE  LEFT  IS  TA  =  TD.  * 


IF  (TA  „  , 

IF  ( TPRME ( I )  .GT.  TC  .AND.  TPRME(I)  .LT.  TA)  THEN 
WPI  =  Z  -  ZPI 
WPF  =  Z  -  L 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN (W2) 

B(I)  =  ROE*N*(BETA**2 .0)*(YY  -  XX) 

END  IF 

IF  (TPRME (I)  .GT.  TD  .AND.  TPRME ( I )  .LT.  TB)  THEN 
WPI  =  Z  -  0.0 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  -  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN (W2) 

B(I )  =  R0E*N*(BETA**2 .0)*(YY  -  XX) 

END  IF 
END  IF 

5  CONTINUE 

END  IF 

x*:k*x;it*****  ***************  ********  A******************************** 

*IN  THE  EVENT  THAT  THE  PREVIOUS  TWO  CONDITIONS  DO  NOT  EXIST,  THIS  * 
*MEANS  WE  ARE  CENTERED  ON  THE  PATH.  HOWEVER,  THE  FOLLOWING  'IF  * 
*STATEMENT  COMPARES  CE  WITH  THETA 1  AND  THETA2  TOGETHER.  THE  WRITE  * 
*STATEMENT  INDICATES  CENTER  ON  THE  SCREEN.  * 

******************************************************************* 

IF  (THETA1  .LT.  CE  .AND,  THETA2  .GE.  CE)  THEN 
WRITE(6 , 113) 

3  FORMAT ( 'PATH  ON  CENTER' ) 

F  =  DTAN(CE) 

ZPC  =  Z  -  (S/F) 

RC  =  SQRT(S**2 .0  +  ((Z  -  ZPC)**2.0)) 
******************************************************************* 
*THE  NEW  VARIABLES  T1  AND  T2  ARE  INTRODUCED.  T1  IS  WHERE  THE  * 

^MINIMUM  OF  THE  U(Z 1 )  CURVE  JUST  TOUCHES  THE  TOP  BOUNDARY  OF  THE  * 
*U1  LINE.  T2  CORRESPONDS  TO  THE  U2  LINE.  T1  IS  WRITTEN  ON  THE  * 
^SCREEN  BY  THE  WRITE  STATEMENT.  * 

******************************************************************* 

T1  =  ((ZPC  +  (BPRME*RC)  -  Ul)/V) 

T2  =  ((ZPC  +  (BPRME*RC)  -  U2)/V) 

WRITER, 333)  Tl/TAU 
3  FORMAT (F10. 4) 

DELR  =  R1  -  R2 
IF  ;  3PP.ME-DELR  .GT.  L)  THEN 
T3  =  TS 
END  IF 

IF  (BPRME*DELR  .LE.  L)  THEN 
T3  =  TD 
END  IF 

TINC  =  (T3  -  Tl)/200 . 0 
DO  30  I  =  1,200 

IF  (I  .EO.  1.0)  THEN 
TPRME (I )  =  T1 
B(I)  =  0.0 
ELSE 

TPRME ( I )  =  T1  +  (REAL(I)*TINC) 


.EO.  TD) 
TPRME (I) 


THEN 


END  IF 

IF  (TPRME(I)  .GE.  T3)  THEN 
B(I )  =  0.0 
GO  TO  100 
END  IF 

A1  =  U1  +  (V*TPRME(I)) 

A2  =  U2  +  (V*TPRME(l)) 

D  =  (BPRME**2.0)*Z  -  A1 
DD  =  (3PRME**2.0)*Z  -  A2 
El  =  ((Z  -  Al}**2.0)  -  51 
E2  =  (( Z  -  A2)**2.0)  -  SI 
IF  (El  .LT.  0.0)  THEN 
El  =  0.0 
END  IF 

IF  (E2  .LT.  0.0)  THEN 
E2  =  0.0 
END  IF 

SRTE1  =  SQRT(El) 

SRTE2  =  SQRT(E2) 

E  =  BPRME*SRTE1 
EE  =  BPRME*SRTE2 

IF  (TPRME(I)  .GT.  T1  .AND.  TPRME(I)  .LT.  T2)  THEN 
ZPM  =  (D  -  E)/Q 
ZPF  =  (D  +  E)/Q 
IF  (ZPH  .LE.  0.0)  THEN 
ZPM  =  0.0 


ZPI  =  ZPM 
END  IF 

IF  (ZPM  .GT.  0.0)  THEN 
ZPI  =  ZPM 
END  IF 

IF  (ZPF  .GE.  L)  THEN 
ZPF  =  L 
END  IF 

IF  (ZPF  .LT.  L)  THEN 
ZPF  =  (D  +  E)/Q 
END  IF 

WPI  *  Z  -  ZPI 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATANW2) 

B(I)  =  ROE*N*(BETA**2.0)*(YY  -  XX) 

END  IF 

IF  (TPRME(I)  .GT.  T2  .AND.  TPRME(I)  .LT.  T3)  THEN 
ZPM  =  (D  -  E)/0 
ZPM2  =  (DD  -  EE)/Q 
IF  (ZPM2  .LE.  0.0)  THEN 


ZPM2  =  (DD  -  EE 
IF  (ZPM2  .LE.  C 


IF  (ZPM2  .LE.  0.0)  THE^ 
ZPF  =  0.0 
END  IF 

IF  (ZPM2  .GT.  0.0)  THEN 
ZPF  =  ZPM2 
END  IF 

IF  (ZPM  .LE.  0.0)  THEN 
ZPI  =  0.0 


END  IF 

IF  'ZPM  .GT.  0.0)  THEN 

7T)T  -  -’t>M 

—  U4  A  1 

END  IF 


WPI  =  Z  -  ZPI 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/S 
YY  =  DATAN(Wl) 

XX  =  DATAN  W2) 

B1  =  ROE*N*(BETA*’I'2.0)*(YY  -  XX) 


ZPI 

ZPF 


(DD  +  EE)/Q 
(D  +  E)/Q  * 
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IF  (ZPI  .GE.  L)  THEN 
ZPI  =  L 
END  IF 

IF  (ZPF  .GE.  L)  THEN 
ZPF  =  L 
END  IF 

IF  (ZPF  .LT.  L)  THEN 
ZPF  =  (D  +  E)/Q 
END  IF 

WPI  =  Z  -  ZPI 
WPF  =  Z  -  ZPF 
W1  =  WPI/S 
W2  =  WPF/ S 
YY  =  DATAN(Wl) 

XX  =  DATAN(W2) 

B2  =  ROE*N*(BETA**2.0)*(YY  -  XX) 

3(1)  =  31  +  B2 
END  IF 

30  CONTINUE 
END  IF 

**  ***  ***  *-k‘k-k'kJckj<jcjcjc7<Jrk  *  *******  **  **********  ***********  ****  ********** 

*ALL  THROUGH  THE  PROGRAM  THE  VALUE  200  FOR  THE  DIMENSIONALIZED  * 
^VARIABLES  TPRME  AND  B  WAS  USED  BECAUSE  IT  BEST  SUITED  THE  ROUTINE* 
*USED  BY  THE  GRAPHICS  PROGRAM  'PLOT1.  THIS  GRAPHICS  ROUTINE  CAN  BE* 
*OBTAINED  FROM  PROFESSOR  GLASS  IN  THE  PHYSICS  DEPARTMENT.  * 

X*XXXXXXX*XXXXX**XXX**XX*XX*XXX*XXXXX*X**X*XX**X*XX**X*X**XX****XXX 

100  DO  50  I  =  1,200 

TPRME ( I )  =  TPRME(I)/TAU 
WRITE (26 , SOO)TPRME(I) ,B(I) 

500  FORMAT (F10 . 5 , F10 . 5) 

50  CONTINUE 


■>vt\ 


APPENDIX  B 
S-Z  PLANE  PROGRAM 


*THIS  PROGRAM  TAKES  THE  DIFFERENT  VARIABLES  AVAILABLE  SUCH  AS  BETA* 
*AND  THE  INDEX  OF  REFRACTION  ETC.  AND  PRODUCES  A  S/L  VS  Z/L  GRAPH  * 
*(S-Z  PLANE).  AGAIN  AS  IN  APPENDIX  A,  THE  GRAPHICS  ROUTINE  USED  * 
*CAN  BE  OBTAINED  BY  PROFESSOR  GLASS  OF  THE  PHYSICS  DEPARTMENT.  * 
******************************************************************* 

DOUBLE  PRECISION  YL ,  YC,  YR,  X,  Z,  BETA,  N,  CO,  RSETME ,  BPRME 
DOUBLE  PRECISION  DELU,  L,  QR,  QL,  QC,  A,  3R,  3L,  BC,  C,  ER,  EL 
DOUBLE  PRECISION  EC,  FR,  FL,  FC,  HR,  HL,  HC,  ZDIVL,  MIN,  INC 
DOUBLE  PRECISION  ONEMIN,  ONEPLS 
INTEGER  I 

DIMENSION  X(200 ) ,  YL(200),  YC(200),  YR(200) 
******************************************************************* 
*THE  VARIABLES  FOR  INPUT  ARE  DESCRIBED  IN  ORDER  TO  BETTER  UNDER-  * 
*STAND  THIS  PROGRAM.  N  IS  THE  INDEX  OF  REFRACTION.  BETA  IS  EQUAL  * 
*TO  V/C  AND  IS  DETERMINED  FROM  THE  LORENTZ  CONTRACTION  FACTOR  * 
*BASED  ON  THE  ENERGY  LEVEL  OF  THE  BEAM.  L  IS  THE  APPROXIMATED  BEAM* 
*LENGTH ,  ZDIVL  IS  THE  MAXIMUM  VALUE  OF  THE  X  AXIS  AND  MIN  IS  THE  * 

*MINIMUM  OF  THE  X  AXIS.  RSETME  IS  THE  TIME  IT  TAKES  FOR  THE  * 

*CURRENT  TO  REACH  THE  MAXIMUM  LEVEL.  * 

*  -k  *  *  -k  *  *  *  *  *  *  *  *  -k  A  *  *  jfc  *  **  *  *  *  *  *  *  -k  "k  •k  "k  *  *  *  *  *  -k  k  k  kk  k  k  k  k  k  k  kk  k  k  k  *  k  k  k  k  k  *  k  *  *  k  k  kk 
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PRINT  *,  'ENTER  N' 

READ  *  N 

PRINT  *,  'ENTER  BETA' 

READ  *  BETA 

PRINT  *,  ‘ENTER  BEAMLENGTH' 

READ  *  L 

PRINT  *,  'ENTER  MAX  ZDIVL' 

READ  *,  ZDIVL 

PRINT  *,  'ENTER  MIN  ZDIVL' 

READ  *  MIN 

PRINT  *,  'ENTER  RISETIME  IN  NSEC' 

READ  *,  RSETME 

INC  =  (ZDIVL  -  MIN)/200 .0 

DO  40  I  =  1,200 

IF  (I  .EQ.  1.0)  THEN 
X(I)  =  MIN 
ELSE 

X(I)  =  MIN  +  REAL(I)*INC 
END  IF 
YR(I)  =  0.0 
YL(I)  =  0.0 
YC(l )  =  0.0 
CONTINUE 
BPRME  =  N*BETA 

DELU  =  RSETME*0 . 29979250*BETA 


ONEMIN  =  1.0  -  (DELU/L) 
ONEPLS  =  1.0  -r  (DELU/L  i 
OR  =  .1.0  -  (DELU/L) ; /BPRME 
QL  =  (1.0  +  (DELU/L) )/ BPRME 
QC  =  1.0 /BPRME 


DO  70  I  «  1,200 
A  =  X(I}**2.0 
BR  =  Q R**2.0 
BL  =  QL**2.0 
BC  =  QC**2 . 0 
C  =  (X(I)  -  1.0)**2.0 
ER  =  (A  -  BR  -  C)**2. 0 


EC  =  (A  -  BC  -  C)**2. 0 
FR  =  ER/(4.0*BR) 

FL  =  EL/(4.0*BL) 

FC  =  EC/(4.0*BC) 

HR  =  FR  -  C 


HL 

=  FL 

-  C 

HC 

=  FC 

-  C 

IF 

(HR  , 

.LT. 

0.0) 

THEN 

HR  = 

:  0.0 

END 

IF 

0.0) 

IF 

ii 

.LT. 

0.0 

THEN 

END 

IF 

0.0) 

IF 

.LT. 

'  n  n 

THEN 

END  IF 

IF  (3PRME  .GT.  ONEMIN)  THEN 
YR(I)  =  SQRT(HR) 

ELSE 

YR(I)  *  0.0 
END  IF 

IF  (BPRME  .GT.  ONEPLS)  THEN 
YL(I)  =  SQRT(HL) 

ELSE 

YL(I)  =  0.0 
END  IF 

IF  (BPRME  .GT.  1.0)  THEN 
YC(I)  =  SQRT(HC) 

ELSE 

.  YC(I)  =  0.0 
END  IF 

70  CONTINUE 

C  *THE  DIMENS IONALIZED  VARIABLES  ARE  USED  AS  FOLLOWS:  X  IS  THE  Z/L  * 
C  *AXIS;  YR  IS  OBTAINED  FROM  SETTING  TB  =  TC;  YL  IS  FROM  TA  =  TD;  * 
C  *YC  IS  FROM  TA  =  TC  OR  TB  =  TD.  * 

Q  ******************************************************************* 

DO  80  I  =  1,200 

WRITE (22, 200)  X(I),  YR(I),  YL(I),  YC(I) 

200  FORMAT (4F10. 5) 

80  CONTINUE 
STOP 
END 
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